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Usually, positions of spacecraft on interplanetary or deep space missions are determined by radar 
tracking from ground stations, a method by which uncertainty increases with distance from Earth. 
As an alternative, a spacecraft equipped with e.g. an X-ray telescope could determine its position 
autonomoulsy via onboard analysis of X-ray pulsar signals. In order to find out which pulsars 
are best suited for this approach and what accuracy can be achieved, we build up a database 
containing the temporal emission characteristics of the ~ 60 X-ray pulsars for which a pulsed 
radiation has been detected by mid 2010. 
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1. Introduction 

The conventional way to obtain data for orbit determination is tracking the spacecraft by 
ground stations on Earth, e.g. with NASA's Deep Space Network (DSN). This technique yields 
very accurate range and range-rate data along the Earth-spacecraft line, but due to limited angular 
resolution large errors can occur in perpendicular directions, resulting in positional errors that grow 
with distance. Interferometric measurements can augment the angular resolution, thereby achiev- 
ing positional accuracies in the order of 4 km per AU of distance between Earth and spacecraft [B]. 
Thus, especially - but not exclusively - for interplanetary and deep space missions, it is desirable 
to have an autonomous navigation solution that works independently from ground stations. 

Pulsars can be used as natural beacons for navigation by comparing pulse arrival times mea- 
sured onboard the spacecraft with predicted arrival times at an inertial reference location - e.g. the 
barycenter of the solar system. The phase difference between the expected and the measured pulse 
arrival corresponds to a run-time difference along the line of sight towards the pulsar and hence 
to a range difference in this direction. Full three-dimensional position information can be deduced 
from the range information along the pulsar lines-of-sight of at least three different pulsars. In 
general, this procedure results in multiple solutions, which can be reduced by either constraining 
the possible spacecraft positions to a finite volume around an initial position assumption, or by 
observing additional pulsars [||, ||] . 



2. Data Analysis and Results 

In order to decide which pulsars are best suited for spacecraft navigation, we re-analyzed all 
pulsar timing data from the X-ray satellites XMM-Newton, Chandra and the ROSSI X-ray Timing 
Explorer. The individual photon arrival times of each data set were corrected for the orbital motion 
of the detector around Earth/Sun via transformation to the solar system barycenter and for the 
orbital motion of the pulsar in case of a binary system. To account for the gradual increase of pulse 
period the spin frequency was modelled as a Taylor expansion, 

f(t) = fo + fo-(t- 1 ) + ifo • (/ - 1 ) 2 , (2. 1) 

with /o, fo, fo being the rotation frequency and its first and second time derivative at some reference 
epoch to, and a pulse number <I> was assigned to each individual arrival time. Since the frequency 



equals the rate of change of pulse number, / = d<J>/df, integration of (2.1) yields 



<&(? ) = $0 + /o • (t - to) + \fo • (t - to) 2 + g/o • (* - to) 3 , (2.2) 

where <I>o denotes the pulse number at to- The pulse numbers are used to generate a mean pulse 
profile reflecting the temporal emission characteristics of the pulsar: (1) For each photon its phase 
<p := (<1> mod 1) is computed, i.e. the fraction of full pulse period at which the photon was detected; 
(2) the domain [0, 1) of is divided into n finite intervals A<f>j := , ^) for i = 1,2, . . . ,n; (3) 
finally, the mean pulse profile is given by the histogram of photon counts per phase interval. The 
optimal value of n depends on the number of recorded photons and the harmonic content of the 
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pulse profile. It can be computed as follows: For each set of arrival times ? ; with i = 1,2, . . . ,N and 
phases 0, := <p(tj) we calculated the statistical variable 

2 / N \ 2-| 

(2.3) 



which is a measure for the periodicity of the signal and evaluated the expression 



//:= max (zl -4m + 4) = Zl -4M + 4 (2.4) 



to get the optimal number M of harmonics [|[]. Then, the formula 

M 



71 = 2.36 



.1/3 

2^ \ m {Zm~Zm-l) 
m=\ 



(2.5) 



with Zq : = yields an estimate for the optimal number of phase intervals Since equation (2.1) 
may not allow to compute f(t) for times r outside the validity range of the ephemerides, the values 
°f fo, fo, fo have to be known at an epoch sufficiently close to the measured photon arrival times. 
We used pulsar ephemerides from the ATNF database ^ for the period folding and produced pulse 
profiles with high signal-to-noise ratios by superposing phase values from several obvservations of 
the same pulsar. 

The navigation algorithm is based on the comparison of pulse arrivals measured onboard the 
spacecraft with those measured at the solar system barycenter. 1 For this purpose, representation of 
the reference profiles by analytical functions is of advantage. Using a least-squares-fit procedure, 
we found approximations of the profiles by either a sum of Gaussians, £,A,exp [— ^ (^^) 2 ] , or of 
sine functions, sin(B,0 +Q) +D, depending on the individual shape of the profile. Concerning 
the use of Gauss or sine functions, the resulting M value of the //-test is a good indicator: For M < 2 
the profile is well approximated by a sum of two or three sine functions, whereas for M > 3 it is 
often better to use a sum of two or more Gaussians. Examples of pulse profiles and their analytical 
representations are shown in Figure [T| These pulse profile templates allow us to measure pulse 
arrivals with high accuracy even for sparse photon statistics - again by using a least-squares fit 
of an adequately adjusted template to the profile in question. The typical error of measuring a 
pulse phase lies in the order of 10~ 3 , which corresponds to an error in range along the pulsar 
line-of-sight of about 1 km for the fastest and 30000 km for the slowest rotating pulsar of our 
database. 2 Evidently, the precision of a pulsar based navigation system strongly depends on the 
choice of pulsars and the accuracy of pulse arrival measurements, wich is subject to the quality of 
the available templates. 

As mentioned above, in order to obtain three-dimensional position information, pulsar timing 
data of at least three different sources have to be collected and analyzed. The spatial arrangement 
of these pulsars is another, significant parameter of the achievable accuracy. The question of how 
pulse arrival errors translate into errors of position in space was addressed as follows: Suppose 



'By pulse arrival we understand the phase of the global maximum in the pulse profile. 

2 PSR B1937+21 with a rotation period P — 1.56 ms and a l-<7 error in phase measurement 8<j> = 2.6 x 1CP 3 ; 
PSR J1841-0456 with P = 11.78 s and 50 = 9.3 x 1(T 3 . The phase error translates into an error of position along the 
pulsar line-of-sight according to the formula Sr = cPS(j), where c is the speed of light. 
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Figure 1: Pulse profiles of PSR B1509-58 (left) and PSR J0633+1746 (right). The data points of the 
pulse profiles can be approximated by a sum of Gauss or sine functions, respectively. The dotted blue lines 
indicate the individual components of the fit (three Gaussians for PSR B1509-58 and two sine functions for 
PSR J0633+1746); their sum is plotted in red. The panels below show the residuals between fit and data 
points, and the red enveloping lines indicate the l-a uncertainties of the binned pulse profile. 



r to be the position of the detector relative to the solar system barycenter (SSB). A pulse arrival 
measured at r will in general be shifted by some value A(p € [0, 1) as compared to a measurement 
at the SSB. Assuming the pulsar signals to be plane waves propagating through flat spacetime, 
this phase shift is given by the formula n-r = (N + A(p)cP, where n is a unit vector from the SSB 
towards the pulsar, N is some integer pulse number, c is the speed of light and P is the pulse period 
at the epoch of measurement. We chose r to be some fixed point in the solar system and for each 
pulsar computed 1) the phase shift associated with that location and 2) a sequence of pulse arrival 
errors, normally distributed around the expected phase shift. The standard deviations of these 
normal distributions are known from our study of pulse arrival measurements. For a combination 
of three different pulsars the last equation can be written in matrix notation and solved for r, 





( nix n\ y n u \ 




r = 


n 2x n 2y n 2z 


l( 




\n 3x n 3y n 3z j 





(N 2 +A((> 2 )cP 2 

\{N 3 + A<$> 3 )cP 3/ 



(2.6) 



Inserting the pulse arrival errors instead of the expected phase shifts in equation ( $J§ gives se- 
quences of position vectors r ( - wich can be compared with the actual position r. A systematic study 



of equation (2.6) for all possbile 3 -combinations from the pulsar set leads to a ranking of pulsar 
combinations. 3 The averaged position errors, i.e. the arithmetic mean of the values |r,- — r\, are 
plotted in Figure || for the first 1000 preferred pulsar combinations. As lower limit, we obtained 
position errors of about 5 km, independent of the location in the solar sytem. However, improve- 
ment in accuracy can be achieved by using pulse profile templates of higher quality. 



The number of 3-combinations from a set of n pulsars is given by the binomial coefficient (") . In our case n = 59 



and hence 



■ 32509 pulsar combinations had to be considered. 
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Figure 2: Position errors for various pulsar combinations. The diagramme shows mean position errors and 
standard deviations for the best 30 combinations (left) and mean position errors for the best 1000 combina- 
tions (right) of three pulsars each. 

3. Summary and Future Prospects 

Currently, we are exploring the achievable accuracy of spacecraft navigation based on X-ray 
pulsar timing data and its feasibility in technical terms. We built up a database of ~ 60 pulsars 
for which a pulsed radiation in X-rays has been detected by mid 2010 (cf. Tables 6.8 and 6.9 in 
Reference [jl|]). The database contains pulsar characteristics like the mean pulse profiles at various 
photon energies, pulsed fractions in various energy bands, pulse profile templates, i.e. approxi- 
mations of pulse profiles by analytical functions and the accuracies of pulse-arrival measurements, 
based on least-squares fitting to the currently available pulse profile templates [j7j]. These data, to- 
gether with the ranking of preferred pulsar combinations, serve as input to the navigation algorithm. 
On the basis of simulated pulse profiles as measured by an arbitrarily moving virtual detector, we 
are studying parameters, such as effective detector area, sensitivity, temporal resolution, etc. As 
to possible implementations, weight, cost, complexity and energy budget of the system, as well as 
technological developments, e.g. light-weighted X-ray mirrors, have to be considered. 
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